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Characterization of uncertainty associated with 
transport-transformation models is often of critical importance, as for 
example in cases where environmental and biological models are employed in 
risk assessment. However, uncertainty analysis using conventional methods 
such as standard Monte Carlo or Latin Hypercube Sampling may not be 
efficient, or even feasible, for complex, computationally demanding models. 

This work introduces a computationally efficient alternative method 
for uncertainty propagation, the Stochastic Response Surface Method (SRSM) . 
The SRSM approximates uncertainties in model outputs through a series 
expansion in normal random variables (polynomial chaos expansion) . The 
unknown coefficients in series expansions are calculated using a limited 

number of model simulations. This method is analogous to approximation of 
a deterministic system by an algebraic response surface. 

Further improvements in the computational efficiency of the SRSM are 
accomplished by coupling the" SRSM with ADIFOR, which facilitates automatic 

calculation of partial derivatives in numerical models coded in Fortran. 
The coupled method, SRSM-ADIFOR, uses the model outputs and their 
derivatives to calculate the unknown coefficients. 

The SRSM and the SRSM-ADIFOR are general methods, and are applicable 
to any model with random inputs. The SRSM has also been implemented as a 
black-box, web-based tool for facilitating its easy use. 

The SRSM and the SRSM-ADIFOR have been applied to a set of 
environmental and biological models. In all the case studies, the SRSM 
required an order of magnitude fewer simulations compared to conventional 
methods, and the SRSM-ADIFOR required even fewer simulations. In addition 
to their computational efficiency, these methods directly provide 
sensitivity information and individual contributions of input uncertainties 
to output uncertainties; conventional methods require substantially larger 
numbers of simulations to provide such information. Thus, the SRSM and the 
SRSM-ADIFOR provide computationally efficient means for uncertainty and 
sensitivity analysis. 

Finally, this research addresses uncertainties associated with model 
.structure and resolution with application to photochemical air quality 
modeling. A three dimensional version of the regulatory Reactive Plume 
Model (RPM), RPM-3D, has been developed and applied to understand model 
uncertainty. 
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The dissertation concerns the extraction, via signal processing, of 
structural information from the scattering of low megahertz, low power 
ultrasonic waves in two specific media of great practical interest--f iber 
reinforced composites and soft biological tissue. 

In fiber reinforced composites, this work represents the first 
measurement of second-order statistics in porous laminates, and the first 
application of Monte Carlo methods to acoustical scattering in composites. 
A numerical model of porous composites backscatter was derived which is 
suitable for direct numerical implementation. The model treats the total 
backscattered field as the result of a two-mode scattering process. In the 
first mode, the void-free composite is treated as a continuously varying 
medium in which the density and compressibility are functions of position. 
The second mode is the distribution of gas voids that failed to escape the 
material before gel, and are dealt with as discrete Rayleigh scatterers. 
Convolution techniques were developed that allowed the numerical model to 
reproduce the long range order seen in the void-free composite. 

The results of the Monte Carlo derivation were coded, and simulations 
run with data sets that duplicate the properties of the composite 
samples used in the study. 

In the area of tissue characterization, two leading methods have been 
proposed to extract structural data from the raw backscattered waveforms. 
Both techniques were developed from an understanding of the periodicities 
created by semi-regularly spaced, coherent scatterers. In the second half 
of the dissertation, a complete analytical and numerical treatment of these 
two techniques was done from a first principles approach. 

Computer simulations were performed to determine the general behavior 
of the algorithms. The main focus is on the envelope correlation spectrum, 
or ECS. Monte Carlo methods were employed, to examine the signal-to-noise 
ratio of the ECS in terms of the variances of the backscattered amplitude, 
phase, spacing, and pulse broadening due to attenuation. Extensive Monte 
Carlo methods have quantified performance limits of the ECS by calculating 

the limits on the S/N ratio in terms of the randomization of 
scattering parameters. A new result is that the S/N ratio was found to be 
controlled by variance of the scattering amplitude and spatial separation 
of the scattering centers. (Abstract shortened by UMI . ) 
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This study is concerned with the classification and segmentation of 
textures. The main emphasis is on the statistical structure of the 
texture rather than its visual structure. The texture is viewed as a 
two-dimensional stochastic spatial process, i.e. a random field. Hence, 
spatial interaction models which capture the spatial structure of random 
fields are used for texture modelling. A class of such models known as 
simultaneous autoregressive (SAR) models is used to represent the texture. 
The estimated parameters of a SAR model fitted to an image are suggested as 
features- for its texture. Results of supervised classification experiments 
using these features indicate that they are powerful. However these 
features are rotation variant, i.e. change if the texture is rotated with 
respect to the camera. 

In order to get a rotation invariant feature set , a new class of 
spatial interaction models called circular autoregressive (CAR) is 
developed. In a CAR model, any pixel value is written as a finite sum of 
intensity values at locations on a circular neighborhood around it and a 
noise sequence. The parameters of this model are rotation invariant. The- 
estimated CAR parameters plus other rotation invariant functions defined on 
SAR parameters are used as textural features. These features have direct 
interpretation in terms of visual attributes of texture. The strong 
discriminating power of the proposed features is shown by using them in 
supervised classification experiments involving differently oriented test 
and training samples from each class. 

A new technique is presented which segments an image' into regions 
of similar texture when no prior information about the textures is 
available. The image is scanned by a small window and a number of textural 
features extracted from each window . ■ Different pairs of features 
corresponding to various windows are plotted. Some of these 2-d plots 
contain distinct clusters corresponding to different texture types. Only 
those two features having the highest clustering power are used. 
Segmentation is done by defining decision boundaries on the clusters 
and mapping them back into the spatial domain. Several experimental results 
are presented. 
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A unique solution is presented for the ancient sampling problem posed 
as follows: One is to test a finite random sample from a large 

population for some property that an unknown fraction of the population 
possesses. Required: To deduce the most accurate possible probabilistic 
statements in the following sense. For a specified probability that the 
unknown fraction be larger than some lower limit , calculate from the 
test result the largest possible lower limit that will have at least this 
probability and the smallest possible upper limit that will have at most 
this probability. Similarly, specify narrowest possible ranges for upper 
limits. These formulas can then be applied to determine for specified lower 
and upper limits two probabilities for their enclosure of the unknown 
fraction, the larger of which is the lowest possible upper limit and the 
smaller of which is the largest possible lower limit. It must be understood 
that these probabilistic statements refer to the conceptual or actual 
results of continued sampling for which the enclosure probabilities are 
held constant and the differing test results give a distribution of lower 
and upper limits. A 'best 1 formula is also given which lies between the 
others and has a zero value of a reasonably defined average error. Use of 
the proposed scheme will reduce needed sample sizes as compared to usual 

statistical procedure in which only 'pessimistic 1 relations are used. 
(Author) 
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large sample example from the Systolic Hypertension in the Elderly 
Program (SHEP) 
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Because many randomized clinical trials study more than one important 
outcome variable, evaluation of efficacy is often difficult and not 
completely satisfactory. This paper considers the use of a procedure for 

endpoint determination described by Follmann et al . , that allows raters 
to integrate subjectively all relevant information about an individual's 
clinical course into a single univariate assessment. To explore the 
method's feasibility, we tested the procedure with data from a completed 
clinical trial, the Systolic Hypertension in the Elderly Program (SHEP) . 
We provided raters blinded to treatment assignment with cards that 
schematically represent the clinical trajectories of SHEP study 
participants. The raters independently ranked these trajectories. The 
method combined ranks across raters to determine a single rank for each 
study participant; we used a rank procedure to test treatment effect. The 
major findings were: (i) the raters showed a high level of concordance of 
rankings; (ii) tests of treatment effect were highly statistically 

significant; (iii) three statistical methods were effective for 
implementing the ranking in the large study size case. These methods were 
use of: (a) scoring rules; (b) incomplete block designs, and (c) 
categorical ranking. 

English Descriptors: Hypertension; Systolic pressure; Human; Elderly; 

Treatment efficiency; Clinical trial; Statistical model; Feasibility; 

Rank statistic ; Statistical test ;. Categorization; Placebo effect 
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French Descriptors: Hypertension arterielle; Pression systolique; Homme; 
Vieillard; Efficacite t-raitement; Essai clinique; Modele statistique; 
Faisabilite; Statistique rang; Test statistique; Categorisation; Effet 
placebo 
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Abstract: For clinical trials with interim analyses, there have been 

methodologies and ■ software to calculate boundaries for comparing 
binomial, normal, and survival data from two treatment groups . 
Jennison & Turnbull (1991) extended Pocock (1977) and O'Brien-Fleming 
(1979) boundaries to t-tests, chi(2)-tests and F-tests for comparing 
normal data from several treatment groups . This paper demonstrates 
that the above boundaries can be applied to a wide variety of test 
statistics based on general parametric settings. We show that 
asymptotically the chi(2) boundaries as well as the corresponding 
nominal significance levels calculated by Jennison & Turnbull can be 
applied to interim analyses based on the score test, the Wald test, and 
the likelihood ratio test for general parametric models. Based on the 
results of this paper, currently available software in group 
sequential testing can be used to calculate the nominal significance 
levels (or boundaries) for group sequential testing based on logistic 
regression, ANOVA, and other parametric methods. 

Descriptors—Author Keywords: interim analysis ; likelihood ratio test ; 
nominal significance level ; partial sum ; score test ; Wald test 
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